Functional renormalization group approach to the sine-Gordon model 
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The renormahzation group flow is presented for the two-dimensional sine-Gordon model within 
the framework of the functional renormalization group method by including the wave-function renor- 
malization constant. The Kosterlitz-Thouless-Berezinski type phase structure is recovered as the 
interpolating scaling law between two competing IR attractive area of the global renormalization 
group flow. 
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/. Introduction. — The two-dimensional (2D) 
Gordon (SG) model, defined by the bare action 



sme- 



S 



[d^^ipf +ucos{P(p) 



(1) 



in Euclidean spacetime, has already received a consider- 
able amount of attention [l|, 0, H, 0) S ISj lZ] since it is 
the simplest non-trivial quantum field theory with com- 
pact variables. This feature is common with non-Abelian 
gauge theories and is supposed to be the key to their con- 
finement mechanism. In two dimensions this is the driv- 
ing force to form a non-trivial phase structure. The SG 
model is known to belong to the universality class of the 
2D Coulomb gas and the 2D-XY spin model which have 
received important applications in condensed matter sys- 
tems, e.g. describing the Kosterlitz-Thouless-Berezinski 
(KTB) [Sj] phase transition of vortices in a thin superfluid 
film. There is a continuous interest in the literature in 
constructing SG type models to understand better the 
vortex dynamics of condensed matter systems [3| ■ 

A more detailed relation between the SG model and 
the XY model in the Villain-approximation is obtained 
by using lattice regularization lO]. The kinetic energy 
is periodic with the same period length as the potential 
energy. Therefore the model supports vortices and has 
a third adjustable parameter, the vortex fugacity z. For 
z — > the vortices are suppressed and the SG model of 
Eq. Il} is recovered in the continuum. The duality trans- 
formation, {P,u,z) — > (27r//3, 2z, u/2) maps the contin- 
uum SG model {z = 0) into the XY model without ex- 
ternal field {u — 0). The Coleman point [Ij, separating 
the renormalizable, asymptotically free phase (/3^ < Stt) 
and the non-renormalizable phase (/3^ > Stt) of the SG 
model is mapped into the KTB point of the XY model. 

The perturbative RG results beyond the local poten- 
tial approximation (LPA) Q can account for the KTB 
phase transition and provide for < Stt in the 

infrared (IR) limit. Recently, by using the flow equation 
approach a different IR limit is obtained for the fre- 
quency, i.e. 47r. However, the latter method is 
not able to recover the leading order perturbative UV re- 
sults for < 4tt, due to the wrong sign of the evolution 



equation derived for the frequency. Functional RG ap- 
proaches have also been used to map the phase structure 
of the SG model but their description is not complete 
since on the one hand the LPA is used 0, 0] and on the 
other hand, the SG model is mapped onto other models 
belonging to the same universality class Therefore, 
the analysis of the SG model is still incomplete. 

Our aim with this work is to determine the complete 
phase structure of the original SG model by extending 
the functional RG analysis beyond the LPA, by includ- 
ing the field-independent wave-function renormalization, 
as well. We use the functional RG method for the ef- 
fective average action [ll|, [13, [3] which enables us to 
treat the wave-function renormalization. The evolution 
arises as the result of the gradual turning on of the field 
fiuctuations according to their increasing amplitude by 
decreasing the control parameter k from the initial value 
A <C fco (with fco the UV cut-off which goes to infinity) 
to zero. 

The phase structure is found to be the global result of 
a competition between an IR fixed line and an IR fixed 
point. The traditional KTB scaling law is actually an 
interpolation between these two effects. 

//. The sine- Gordon model. — The functional renor- 
malization group equation for the effective action of an 
Euclidean field theory is [lH 



, ^ ^ 1^ kdkRk 



(2) 



where the notation ' = d/dip is used and the trace Tr 
stands for the integration over all momenta. We use a 
power-law type regulator function 



Rk 



P 



(3) 



with the parameter b > 1. Eq. ([2]) has been solved over 
the functional subspace defined by the ansatz 



Tk = 



Vk{iPx) 



(4) 
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with the local potential Vk{(p) = J2'^=i Un{k) cos{n(p) 
and the field-independent wave-function renormalization 
z{k). Eq. ^ leads to the evolution equations [3] 



dkVk 



1 



VkkdkRk, 



kdkz = VoVl"^ [ VlkdkRi 
Jp 



(5) 



with Vk = + R^ + VI') and Vo = (2^)-i J^^ 

being the projection onto the field-independent subspace. 

///. Linearized scaling at the Coleman point. — We as- 
sume A^^k^^\V^'\, keep the leading order terms in V^' 
in the Taylor-expansion of the r.h.s of Eqs. ([5]) and ([6]) 
and retain a single Fourier mode in the potential Vk 
for simplicity. The two-dimensional momentum integrals 
can easily be performed, giving 



{2 + kdk)ui = - — ui, 
47rz 

kdkZ = — - 



,2-2/b 



(7) 
(8) 



where the dimensionless couplings Un — k '^Un are intro- 
duced, and 



(9) 



The UV evolution Eqs. Q and ^ clearly show that the 
critical value z* = I/Stt at the Coleman point is inde- 
pendent of the blocking parameter b. Furthermore the 
sharp cutoff limit b oo gives infinite value in the r.h.s. 
of Eq. ([8]) signalling the impossibility of introducing the 
wave-function renormalization in that case. The RG tra- 
jectories obtained by integrating Eqs. ^ and 



uiiz) 



(8^)l-2/bcb 
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indicate turning points in the vicinity of the fixed point, 
at {u\ = £ti(z*), z*). Such a flow exhibits the well-known 
features of the KTB type phase transition. Actually we 
see the dual of that transition as explained in the Intro- 
duction. 

Thus Eqs. ([7]) and ^ provide similar evolution around 
the KTB fixed point as the one already obtained by a 
perturbative RG analysis Q and the flow equation ap- 
proach [4| for the SG model, and also by the real-space 
RG for the two-dimensional Coulomb gas . The KTB 
phase transition is characterized by the exponential de- 
pendence of the correlation length on the inverse of the 
square-root of the reduced temperature t cx . The cor- 
relation length ^ can be read off from the scale k* ^ 1/^ 
where the RG trajectories show up their turning points. 
Inserting back the solution PH)) into Eq. ([5]) one obtains 

^ gV^/(fil8v/S?)+S-(f,-l)(2&-l)v^2i-'3/'',75/=^-2/''6-2 + Ci(fij") 

(11) 



which is the typical scaling law for KTB type phase tran- 
sitions, modified by analytic corrections vanishing for 
u\ —t 0. It is worthwhile mentioning that only the quan- 
titative details depend on the choice of the parameter b 
in the formula PT|) . Using Eqs. ([7]) and ^ the critical 
exponent can also be calculated via the vortex-vortex 
correlation function Ts*] and it is proved to take the value 
?7 = 1/4 independently of the parameter h. 

IV. Coleman point, revisited as the dual KTB point. — 
Let us now take into account the higher-order terms of 
the Taylor expansion in V^' , as well as the higher har- 
monics of the local potential 0]. We choose 6=1, corre- 
sponding to the Callan-Symanzik RG scheme fisj which 
is free of UV divergences for d = 2 and ultralocal. The 
evolution equations assume a simpler form rendering eas- 
ier the handling of the higher Fourier modes. 

The Fourier transform of Eqs. (O and ([6]) produces a 
set of coupled equations for £t„ Q and z. We refer to the 
solution of these equations with 10 Fourier modes as the 
full solution. According to our experience the retaining of 
more Fourier harmonics modifies the flow in a negligible 
manner. By restricting the solution to a single Fourier 
mode one obtains the evolution equations 



{2 + kdk)ui - 
kdkZ = 



1 



2ttuiz 
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(12) 
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whose solution will be referred as exact in Hi. The first 
two terms in the Taylor expansion of Eq. (|12p can be 
identified with the approximation used in ^] by a proper 
transformation of the parameters in Callan-Symanzik RG 
scheme. The higher-order terms make negligible effect on 
the evolution even in the neighbourhood of the turning 
point {ul, z*) of the RG trajectory. We refer to the solu- 
tion of Eqs. ^ and ([8]) as the Zmearized solution. The RG 
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FIG. 1: Phase diagram of the SG model. The solid/ 
dashed/dotted lines show the RG trajectories for the lin- 
earized/exact in Ml /full solutions, respectively. The wide solid 
line depicts the separatrix. 



trajectories are plotted in Fig.[Tl they move to the left as 
k is decreased. This picture is reminiscent of the usual 
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KTB phase structure. What we see here is actually the 
vicinity of the dual of the KTB point of the XY model 
[lot . One can see that the higher harmonics modify the 
RG trajectories rather slightly as compared to both the 
linearized and the exact in ui solutions. Although the 
Coleman point lies at a crossover scale between the UV 
and IR scaling regions for the trajectories above the sep- 
aratrix, the UV scaling remains valid in its vicinity. 

Note the smallness of the z-interval covered. The z- 
dependence of Eq. ^ is weak under the two solid lines of 
the separatrix, boardering three different regions, where 
the regions from the left to the right correspond to the 
renormalizable, the non-renormalizable and the asymp- 
totically free regimes of the SG model, respectively. 

V. Non-perturbative scaling at the crossover. — Let us 
now turn our attention to the new phase which is opened 
up by the evolution of the wave function renormalization 
constant in the middle of the figure, above the separatrix. 
The z-dependence is crucial here, it prevents the system 
to come to standstill where the renormalized trajectory 
is stationary in ui, just above the Coleman point. It 
has been established that all Fourier modes are irrelevant 
(decrease with k) before the crossover and they turn to 
relevant (starts to increase as k is further decreased) at 
A;*, showing very weak n-dependence at the location of 
the turning point, fc* « k* . Wefindz(A:*) = z* = l/(87r), 
in a manner similar to the case z — 1 [7]. 

The vertical line z = z* appears to be a single IR 
stable fixed point as far as the evolution of the poten- 
tial is considered only. In fact, the values of the coupling 
constants, u* = Un(k*), determined by the Fourier trans- 
form of the evolution equation Eq. ^ satisfy at this line 
the condition that the ratios 



(14) 



~,*2n ' 



are universal constants, C2 = 1/12, C3 = 1/96, C4 = 
13/8640, C5 = 97/414720, etc. Q. We recover renor- 
malizability and asymptotical freedom in this phase be- 
cause the dynamics is characterized by a single coupling 
strength, ui > at and below the crossover scale. This 
is a non-perturbative phenomenon because the crossover 
"fixed point" is not Gaussian. 

We see furthermore the subtle meaning of the " KTB 
fixed point". As soon as one goes beyond the LPA the 
Coleman point ceases to be a fixed point and is separating 
different phases only under the separatrix of Fig.[Tl where 
the beta functions have a common analytic structure Q 
and something irregularity shows up in the deep IR region 
of the symmetry broken phase only. 

A more detailed and explicit similarity with the KTB 
scaling of the XY model is found by introducing the cor- 
relation length ^ by identifying it with the inverse cut- 
off at the crossover. The numerical results are shown 
in Fig. [51 The various approximations, i.e. the full 
solution, the solution exact in ui, and the linearized 



one give the same critical behaviour as Eq. (fTTj) . show- 
ing that neither the inclusion of higher-order terms in 
Eqs. ^ and ([8|) nor that of the higher harmonics affect 
the type of the phase transition, their effects are negli- 
gible. The reduced temperature is given formally by the 
wave-function renormalization constant at the UV cutoff 
as t = [z(A) - Zs{A)]/zs{A), where (l/8TrZs{A),ul{A)) is 
a point of the separatrix. The turning point u* is shown 
in the inset of Fig. [2] as the function of the reduced tem- 
perature t for the linearized solution, giving 



uf = qt + 0{t^) 



(15) 



as in the XY model [8,], and in the Coulomb gas [l5|. 
The same relation is recovered for the exact solution in 
ill and for the full solution, as well. The critical scaling 
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FIG. 2: The inverse of the turning value ul of the fundamental 
amplitude is plotted versus the correlation length ^ ~ 1/fc* 
for the linearized, exact in ui and full solutions, denoted by 
different point types. The slope of the solid line is -y/S/STr^ ~ 
0.52 according to Eq. (Illf) . In the inset the turning value ul 
is plotted against the reduced temperature t. The slope of 
the fitted dashed line is 0.5. 



relations ([TT|) and signal directly that there is a KTB 
type phase structure in the SG model. 

VI. The IR scaling regime. — The trajectories end in a 
line of Gaussian IR fixed point in the non-renormalizable 
phase. All coupling strengths of the potential are irrel- 
evant and this implies that the evolution of the wave- 
function renormalization z is extremely weak. The LPA 
can be used and the well-known IR scaling is recovered, 
including the unusual feature of the non-availability of 
the concept of relevant or irrelevant operators . There 
is a line of Gaussian fixed points in the asymptotically 
free phase, too, but these fixed points are UV and their 
scaling laws are linearizable. The IR scaling is diffi- 
cult to establish numeri cally because of the instability 
of the Fourier expansion [l^ in any RG scheme, used so 
far. Nevertheless it is unambiguous from numerics that 
z tends to be big as the scale k is decreased, while ui 
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FIG. 3; The phase structure of the SG model given by the 
numerical solution of Eqs. (|12|l and (|13|l . 



remains finite. One can get a clear picture of the phase 
structure of the IR scaling regime by omitting the effect 
of the higher harmonics, a frequently used approxima- 
tion 0, After introducing to = Vl — u^, X = '^1 zu) 
and dt — u'^kdk we arrive at the evolution equations 

dtX - X^i^-2x(l--^) + ^(l-c.), (16) 

possessing two lines of Gaussian fixed points separated 
by the well-known Coleman (alias KTB) fixed point, 
(l/87rz,{ti) = (1,0), and an additional (IR) fixed point 
(l/87rz,{ti) — (0, 1) (see Fig.[3]for the RG trajectories). 

Such a modification of the scaling laws which is be- 
lieved to preserve the qualitative features of the RG flow 
makes the IR fixed points explicit in the complete phase 
diagram. It also demonstrates that the hyperbolic na- 
ture of the flow in the vicinity of the KTB-Coleman point 
stems from global effects, the competition between two 
regions of the phase diagram. This is the attraction of 
the line of Gaussian IR fixed points of the symmetrical 
phase, dominated by the kinetic energy on the one hand 
and of the non-Gaussian, IR fixed point which is domi- 
nated by the potential energy on the other. 

The feature lost in this approximation is that the ef- 
fective potential, built up in the IR region makes the 
evolution equation singular, which is typical for phases 
with spontaneously broken symmetry [17| , automatically 
guaranteeing the superuniversality for the potentials. In 
fact, the precise treatment with no expansion should give 
the superuniversal potential 14^0 — ~\4>'^ due to the 



Maxwell cut QHIli and l/z{k ^ 0) = 0. 
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